Three years of experience with the STELLA robotic observatory

Since May 2006, the two STELLA robotic telescopes at the Izana observatory in Tenerife, Spain, delivered an almost uninterrupted stream of scientific data. To achieve such a high level of autonomous operation, the replacement of all troubleshooting skills of a regular observer in software was required. Care must be taken on error handling issues and on robustness of the algorithms used. In the current paper, we summarize the approaches we followed in the STELLA observatory

3D simulations of rising magnetic flux tubes in a compressible rotating interior: The effect of magnetic tension

Context: Long-term variability in solar cycles represents a challenging constraint for theoretical models. Mean-field Babcock-Leighton dynamos that consider non-instantaneous rising flux tubes have been shown to exhibit long-term variability in their magnetic cycle. However a relation that parameterizes the rise-time of non-axisymmetric magnetic flux tubes in terms of stellar parameters is still missing. Aims: We aim to find a general parameterization of the rise-time of magnetic flux tubes for solar-like stars. Methods: By considering the influence of magnetic tension on the rise of non-axisymmetric flux tubes, we predict the existence of a control parameter referred as $\Gamma_{\alpha_1}^{\alpha_2}$. This parameter is a measure of the balance between rotational effects and magnetic effects (buoyancy and tension) acting on the magnetic flux tube. We carry out two series of numerical experiments (one for axisymmetric rise and one for non-axisymmetric rise) and demonstrate that $\Gamma_{\alpha_1}^{\alpha_2}$ indeed controls the rise-time of magnetic flux tubes. Results: We find that the rise-time follows a power law of $\Gamma_{\alpha_1}^{\alpha_2}$ with an exponent that depends on the azimuthal wavenumber of the magnetic flux loop. Conclusions: Compressibility does not impact the rise of magnetic flux tubes, while non-axisymmetry does. In the case of non-axisymmetric rise, the tension force modifies the force balance acting on the magnetic flux tube. We identified the three independent parameters required to predict the rise-time of magnetic flux tubes, that is, the stellar rotation rate, the magnetic flux density of the flux tube, and its azimuthal wavenumber. We combined these into one single relation that is valid for any solar-like star. We suggest using this generalized relation to constrain the rise-time of magnetic flux tubes in Babcock-Leighton dynamo models.Comment: 18 pages, 15 figures, 6 tabula

On the binary orbit of Henry Draper one (HD 1)

We present our final orbit for the late-type spectroscopic binary HD 1. Employed are 553 spectra from 13 years of observations with our robotic STELLA facility and its high-resolution echelle spectrograph SES. Its long-term radial-velocity stability is $\approx$50m/s . A single radial velocity of HD 1 reached a rms residual of 63m/s, close to the expected precision. Spectral lines of HD 1 are rotationally broadened with $v\sin i$ of 9.1$\pm$0.1 km/s . The overall spectrum appears single-lined and yielded an orbit with an eccentricity of 0.5056$\pm$0.0005 and a semi-amplitude of 4.44km/s . We constrain and refine the orbital period based on the SES data alone to 2318.70$\pm$0.32d, compared to 2317.8$\pm$1.1d when including the older data set published by DAO and Cambridge/Coravel. Owing to the higher precision of the SES data, we base the orbit calculation only on the STELLA/SES velocities in order not to degrade its solution. We redetermine astrophysical parameters for HD 1 from spectrum synthesis and, together with the new {\it Gaia} DR-2 parallax, suggest a higher luminosity than published previously. We conclude that HD 1 is a slightly metal-deficient K0 III-II giant 217 times more luminous than the Sun. The secondary remains invisible at optical wavelengths. We present evidence for the existence of a third component

On the lithium abundance of the visual binary components ξ Boo A (G8V) and ξ Boo B (K5V)

A spectroscopic investigation of the lithium resonance doublet in ξ Boo A and ξ Boo B in terms of both abundance and isotopic ratio is presented. We obtained new R = 130,000 spectra with a signal-to-noise ratio (S/N) per pixel of up to 3200 using the 11.8 m LBT and PEPSI. From fits with synthetic line profiles based on 1D-LTE MARCS model atmospheres and 3D-NLTE corrections, we determine the abundances of both isotopes. For ξ Boo A, we find A(Li) = 2.40 ± 0.03 dex and 6Li/7Li <1.5 ± 1.0% in 1D-LTE, which increases to ≈2.45 for the 3D-NLTE case. For ξ Boo B we obtain A(Li) = 0.37 ± 0.09 dex in 1D-LTE with an unspecified 6Li/7Li level. Therefore, no 6Li is seen on any of the two stars. We consider a spot model for the Li fit for ξ Boo B and find A(Li) = 0.45 ± 0.09 dex. The 7Li abundance is 23 times higher for ξ Boo A than the Sun's, but three times lower than the Sun's for ξ Boo B while both fit the trend of single stars in the similar-aged M35 open cluster. Effective temperatures are redetermined from the TiO band head strength. We note that the best-fit global metallicities are −0.13 ± 0.01 dex for ξ Boo A but +0.13 ± 0.02 dex for ξ Boo B. Lithium abundance for the K5V benchmark star 61 Cyg A was obtained to A(Li) ≈ 0.53 dex when including a spot model but to ≈0.15 dex without a spot model

Flares in Open Clusters with K2. I. M45 (Pleiades), M44 (Praesepe) and M67

The presence and strength of a stellar magnetic field and activity is rooted in a star's fundamental parameters such as mass and age. Can flares serve as an accurate stellar "clock"? To explore if we can quantify an activity-age relation in the form of a flaring-age relation, we measured trends in the flaring rates and energies for stars with different masses and ages. We investigated the time-domain photometry provided by Kepler's follow-up mission K2 and searched for flares in three solar metallicity open clusters with well-known ages, M45 (0.125 Gyr), M44 (0.63 Gyr), and M67 (4.3 Gyr). We updated and employed the automated flare finding and analysis pipeline Appaloosa, originally designed for Kepler. We introduced a synthetic flare injection and recovery subroutine to ascribe detection and energy recovery rates for flares in a broad energy range for each light curve. We collected a sample of 1 761 stars, mostly late-K to mid-M dwarfs and found 751 flare candidates with energies ranging from $4\cdot10^{32}$ erg to $6\cdot10^{34}$ erg, of which 596 belong to M45, 155 to M44, and none to M67. We find that flaring activity depends both on $T_\mathrm{eff}$, and age. But all flare frequency distributions have similar slopes with $\alpha \approx2.0-2.4$, supporting a universal flare generation process. We discuss implications for the physical conditions under which flares occur, and how the sample's metallicity and multiplicity affect our results.Comment: 17 pages, 11 figures, appendix. Accepted to A&

Stellar rotation, binarity, and lithium in the open cluster IC4756

An important aspect in the evolutionary scenario of cool stars is their rotation and the rotationally induced magnetic activity and interior mixing. Stars in open clusters are particularly useful tracers for these aspects because of their known ages. We aim to characterize the open cluster IC4756 and measure stellar rotation periods and surface differential rotation for a sample of its member stars. Thirty-seven cluster stars were observed continuously with the CoRoT satellite for 78 days in 2010. Follow-up high-resolution spectroscopy of the CoRoT targets and deep Str\"omgren $uvby\beta$ and H$\alpha$ photometry of the entire cluster were obtained with our robotic STELLA facility and its echelle spectrograph and wide-field imager, respectively. We determined high-precision photometric periods for 27 of the 37 CoRoT targets and found values between 0.155 and 11.4 days. Twenty of these are rotation periods. Twelve targets are spectroscopic binaries of which 11 were previously unknown; orbits are given for six of them. Six targets were found that show evidence of differential rotation with $\Delta\Omega/\Omega$ in the range 0.04-0.15. Five stars are non-radially pulsating stars with fundamental periods of below 1d, two stars are semi-contact binaries, and one target is a micro-flaring star that also shows rotational modulation. Nine stars in total were not considered members because of much redder color(s) and deviant radial velocities with respect to the cluster mean. H$\alpha$ photometry indicates that the cluster ensemble does not contain magnetically over-active stars. The cluster average metallicity is -0.08$\pm$0.06 (rms) and its logarithmic lithium abundance for 12 G-dwarf stars is 2.39$\pm$0.17 (rms). [...]Comment: A&A, in pres

Time-series photometric SPOT modeling. II. Fifteen years of photometry of the bright RS CVN binary HR 7275.

We present a time-dependent spot modeling analysis of 15 consecutive years of V-band photometry of the long-period (Porb = 28.6 days) RS CVn binary HR 7275. This baseline in time is one of the longest, uninterrupted intervals a spotted star has been observed. The spot modeling analysis yields a total of 20 different spots throughout the time span of our observations. The distribution of the observed spot migration rates is consistent with solar-type differential rotation and suggests a lower limit of the differential-rotation coefficient of 0.022 +/-0.004. The observed, maximum lifetime of a single spot (or spot group) is 4.5 years, the minimum lifetime is approximately one year, but an average spot lives for 2.2 years. If we assume that the mechanical shear by differential rotation sets the upper limit to the spot lifetime, the observed maximum lifetime in turn sets an upper limit to the differential-rotation coefficient, namely 0.04 +/- 0.01. This would be differential rotation just 5 to 8 times less than the solar value and one of the strongest among active binaries. We found no conclusive evidence for the existence of a periodic phenomenon that could be attributed to a stellar magnetic cycle